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r~| ! Abstract. A dressed-quark core contribution to nucleon electromagnetic form 

"V* ' factors is calculated. It is defined by the solution of a Poincare covariant Fad- 

Q ■ deev equation in which dressed-quarks provide the elementary degree of free- 

^ ! dom and correlations between them are expressed via diquarks. The nucleon- 

I ^ i ' photon vertex involves a single parameter; i.e., a diquark charge radius. It is 

argued to be commensurate with the pion's charge radius. A comprehensive 

, analysis and explanation of the form factors is built upon this foundation. A 

! particular feature of the study is a separation of form factor contributions into 

' those from different diagram types and correlation sectors, and subsequently 

. a flavour separation for each of these. Amongst the extensive body of results 

^ I that one could highlight are: r"'" > r"'"^, owing to the presence of axial-vector 

T— I ' quark-quark correlations; and for both the neutron and proton the ratio of 

00 . Sachs electric and magnetic form factors possesses a zero. 
O ■ 

> ' 

^ ' 1 Introduction 



Owing in part to the relatively simple nature of the virtual photon as a probe, a 
reliable explanation of electromagnetic form factors provides information on the 
distribution of a nucleon's characterising properties; e.g., total- and angular- 
momentum, amongst its QCD constituents. Since contemporary experiments 
employ > M^; i.e., momentum transfers in excess of the nucleon's mass, a 
veracious understanding of the body of extant data requires a Poincare covariant 
description of the nucleon. Poincare covariance and the vector exchange nature 
of QCD guarantee the existence of nonzero quark orbital angular momentum in 
a hadron's rest-frame bound-state amplitude [HE]. 

In fact the challenge is compounded owing, e.g., to the running of the dressed- 
quark mass O m m El El E] . This entails that a quantum field theoretic treatment 
of hadron structure and electromagnetic interactions is generally necessary in 
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order to provide understanding in terms of QCD's genuine elementary degrees of 
freedom. The dressed light-quark mass function at infrared momenta is roughly 
100-times larger than the current-quark mass. This marked enhancement is a 
corollary of dynamical chiral symmetry breaking (DCSB) and owes primarily to 
a dense cloud of gluons that clothes a low- momentum quark [9]. (The dressing 
gluons also acquire mass dynamically [lO].) It means that the Higgs mechanism is 
largely irrelevant to the bulk of normal matter in the universe. Instead the single 
most important mass generating mechanism for light-quark hadrons is the strong 
interaction effect of DCSB; e.g., one can identify DCSB as being responsible for 
98% of a proton's mass. It has long been argued that form factors are a sensitive 
probe of this effect [H] . 

Recent years have seen rapid experimental and theoretical progress in the 
study of nucleon electromagnetic form factors, which is reviewed, e.g., in Refs. \12\ 
[13]. Despite this, questions remain unanswered, amongst them: can one formulate 
an impulse-like approximation for hadron form factors and, if so, in terms of 
which degrees of freedom; is there a valid mapping of form factors into statements 
about the distribution of charge and magnetisation within a nucleon; and what 
role is played by pseudoscalar mesons in hadron electromagnetic structure and 
can one describe this in a quantitative, model-independent fashion? Herein we 
contribute to the discussion of these issues. 

In Sect. [2] we recapitulate briefly upon a Poincare covariant Faddeev equa- 
tion for the nucleon, in which the primary element is the dressed-quark with 
its strongly momentum dependent mass function. The Faddeev equation solu- 
tion defines a nucleon's dressed-quark core. The study of baryons in this way 
sits squarely within the ambit of the application of Dyson-Schwinger equa- 
tions (DSEs) in QCD [T3]. Since the DSEs admit a nonperturbative symmetry- 
preserving truncation scheme \15\ [T6l \T7\ I18j . which e.g. has enabled the proof 
of numerous exact results for pseudoscalar mesons [19l [20l [211 [2l] i the approach 
holds particular promise as a means of unifying the treatment of meson and 
baryon observables that preserves all global and local corollaries of DCSB with- 
out fine-tuning [23]. The coupling of a photon to the nucleon's dressed-quark 
core is detailed in Sect.[3| 

In Sect. [4] we discuss the interpretation of form factors and present a perspec- 
tive on the circumstances under which the three dimensional Fourier transform 
of a Breit-frame Sachs form factor can reasonably be understood in terms of a 
charge or magnetisation density. 

Sections [5HZI are extensive. They detail our computed results and the un- 
derstanding they provide. All electromagnetic form factors of the proton and 
neutron are described along with their decomposition into individual flavour, di- 
agram and diquark contributions, the meaning of which will subsequently become 
apparent. 

We consider form factor contributions arising from pseudoscalar meson loops 
in Sect. [8] and exemplify the manner in which they add to the dressed-quark core 
results. We wrap-up in Sect M 
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Figure 1. Poincare covariant Faddeev equation, Eq. HA.11[) . employed herein to calculate nu- 
cleon properties. ^ in Eq. (|A.l|l is the Faddeev amplitude for a nucleon of total momentum 
P = Pq + Pd- It expresses the relative momentum correlation between the dressed-quark and 
-diquarks within the nucleon. The shaded region demarcates the kernel of the Faddeev equa- 
tion, Sect. lA.2l in which; the single line denotes the dressed-quark propagator, Sect. lA.2!T] F 
is the diquark Bethe-Salpeter-like amplitude, Sect. lA.2!2l and the double Ime is the diquark 
propagator, Sect. lA.2!3l 



2 Nucleon Model 



In quantum field theory a nucleon appears as a pole in a six-point quark Green 
function. The pole's residue is proportional to the nucleon's Faddeev amplitude, 
which is obtained from a Poincare covariant Faddeev equation that adds-up all 
possible quantum field theoretical exchanges and interactions that can take place 
between three dressed-quarks. Canonical normalisation of the Faddeev amplitude 
guarantees unit residue for the s-channel nucleon pole in the = three- 
quark vacuum polarisation diagram and entails unit charge for the proton. 

A tractable truncation of the Faddeev equation is based |24j on the obser- 
vation that an interaction which describes mesons also generates diquark corre- 
lations in the colour-3 channel [25]. The dominant correlations for ground state 
octet and decuplet baryons are scalar (O"*") and axial- vector (1"*") diquarks be- 
cause, for example, the associated mass-scales are smaller than the baryons' 
masses [26l [27j , namely (in GeV) 



m 



U(I\ 



0+ 



0.7-0.8, m 



(uu) 



1+ 



m 



m 



0.9- 1.0. 



The kernel of the Faddeev equation is completed by specifying that the quarks 
are dressed, with two of the three dressed-quarks correlated always as a colour-3 
diquark. As illustrated in Fig[Tl binding is then effected by the iterated exchange 
of roles between the bystander and diquark-participant quarks. 



The Faddeev equation that we employ is explained in Appendix A: Faddeev 
Equation. With all its elements specified, as described therein, the equation can 
be solved to obtain the nucleon's mass and amplitude. Owing to Eq. (|A.34p . in 
this calculation the masses of the scalar and axial-vector diquarks are the only 
variable parameters. The axial-vector mass is chosen so as to obtain a desired 
for the Z\0 and the scalar mass is subsequently set by requiring a particular 



mass : 



^This is natural because the spin- and isospin-3/2 A contains only an axial-vector diquark. The 
relevant Faddeev equation is not different in principle to that for the nucleon. It is described 
in Ref. [30]. 
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nucleoli mass. 

We have written here of desired rather than experimental mass values because 
it is known that the masses of the nucleon and A are materially reduced by 
pseudoscalar meson loop effects. This is detailed in Refs. [281 [29]. Hence, a baryon 
represented by the Faddeev equation described above must possess a mass that 
is inflated with respect to experiment so as to allow for an additional attractive 
contribution from the pseudoscalar mesons. As in previous work [301 [3T| [32| [33] 
and reported in Table [H we require Mn = 1.18 GeV and = 1.33 GeV. 
The results and conclusions of our study are essentially unchanged should even 
larger masses and a smaller splitting — M^r be more realistic, a possibility 
suggested by Refs. j23[ 134] . The relationship between the A-N mass splitting 
and that between the axial-vector and scalar diquark correlations is sketched in 
Ref. [35]. 



3 Nucleon Electromagnetic Current 

The nucleon's electromagnetic current is 

MP',P) = ieu{P')A^{q,P)u{P), (2) 
= ien(P') (7M^i(Q') + ^cT^.Q.i"2(Q'))u(P), (3) 

where P {P') is the momentum of the incoming (outgoing) nucleon, Q = P' — P, 
and Fi and F2 are, respectively, the Dirac and Pauli form factors. They are the 
primary calculated quantities, from which one obtains the nucleon's electric and 
magnetic (Sachs) form factors 

Ge{Q^) = Fi{Q^) - ^F2{Q') , Gm{Q') = FiiQ^) + F2{Q^) . (4) 

Static electromagnetic properties are associated with the behaviour of these 
form factors in the neighbourhood of ~ 0. The nucleons' magnetic moments 
are defined through 

^Jin = ^^n = , ^Jip = I + Kp = GljiS)) , (5) 



Table 1. Mass-scale parameters (in GeV) for the scalar and axial- vector diquark correlations, 
fixed by fitting nucleon and A masses offset to allow for "pion cloud" contributions [28] ■ We also 
list ujjp = —T^rujp, the width-parameter in the {qq)jp Bethe-Salpeter amplitude, Eqs. (|A.28|I 



& ([A:29| : its inverse is an indication of the diquark's matter radius. Row 3 illustrates effects of 
omitting the l+-diquark correlation: the A cannot be formed and Mm is significantly increased. 
Evidently, the 1^-diquark provides significant attraction in the Faddeev equation's kernel. 
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0.56=1/(0.35 fm) 
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Figure 2. Vertex which ensures a conserved current for on-shell nucleons described by the 
Faddeev amphtudes, ^ij, described in Sect .[2] and [Appendix A:] Faddeev Equation. The single 
hue represents S(p), the dressed-quark propagator, Sec. lA.2."n and the double line, the diquark 
propagator, Sec. lA.2.^ F is the diquark Bethe-Salpeter amplitude. Sec. I A. 2.^ and the remaining 
vertices are described in |Appendix C:| the top-left image is Diagram 1; the top-right. Diagram 2; 
and so on, with the bottom-right image. Diagram 6. 



where kn, N = n,p, are referred to as the anomalous magnetic moments; and 
the electric and magnetic rms radii via 



s=0 



n) 



d 



-6-lnGj^(s) 



(6) 
(7) 



In order to calculate the electromagnetic form factors one must know the 
manner in which the nucleon described in Sect. [2] couples to a photon. That is 



derived in Ref. [36], illustrated in Fig. [2] and detailed in Appendix C: Nucleon- 
Photon Vertex. As apparent in that Appendix, the current depends on the elec- 
tromagnetic properties of the diquark correlations. 

Estimates exist of the size of diquark correlations. For example, a first Fad- 
deev equation study of nucleon form factors [37] found a scalar diquark radius of 



'Mo4 



0.8 Ttt, where r^r is the pion charge radius within the same model. One 



obtains a similar result in a DSE calculation [38] that provides a good description 
of pseudoscalar and vector meson properties; i.e., 



0.7 fm. 



''{ud) 



1+ 



0.8 fm. 



(8) 
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where the last result is an esthnate based on the p-meson/vr-meson radius-ratio 
[39\ I40j . From another perspective, numerical simulations of quenched lattice- 
regularised QCD suggest a scalar-diquark matter-radius 



rf,i =l.l±0.2fm. (9) 

Mo+ ^ ' 

It is thus evident that diquark correlations within a baryon are not pointlike. 
Hence, with increasing ^ interaction diagrams in which the photon resolves a 
diquark's substructure must be suppressed with respect to contributions from di- 
agrams that describe a photon interacting with a bystander or exchanged quark. 
These latter are the only hard interactions with dressed-quarks allowed in a nu- 
cleon. One can therefore improve on Refs. |3H I32j by introducing a diquark form 
factor. This is expressed in Eqs. (fCTTS]) . (fOlil) and (fOMj) . 

We use a one-parameter dipole because the system involves two quarks. The 
parameter is a length-scale that characterises the diquark radius. In the absence 
of an explicit calculation of the axial- vector diquark's radius, we employ the 
same value for scalar and axial- vector diquarks. Owing to differences between the 
formulation of our nucleon model and the DSE truncation employed in Ref. [38], 
the values quoted in Eq. ([8]) provide only a loose constraint on this parameter. 
It's value does not have a large effect on form factors for < 2GeV^ but does 
influence their evolution thereafter. For example, it influences the position of the 
zero in G^{Q'^): a larger diquark radius shifting the zero further from the origin. 
Computations have been analysed with four values: r^g = 0.0, 0.4, 0.8 and 1.2 fm. 
Unless otherwise stated, the results reported herein were obtained with 

rgg = 0.8fm. (10) 



4 On Interpreting Form Factors 

Now that the Faddeev equation and a consistent Ward-Takahashi-identity con- 
serving current are completely defined, the calculation of a nucleons' electromag- 
netic form factors is a straightforward numerical exercise. However, in light of 
Refs. |42| I43j we judge it worthwhile to comment on their putative interpretation 
in terms of charge and magnetisation densities before presenting our results. 

Such an interpretation rests on the existence of a quantitatively reliable ex- 
pression for the form factors in terms of a current in which the interacting con- 
stituents are well-defined and distinct, for then the charge and current carrying 
quanta are unambiguous. This is achieved through a current of impulse approxi- 
mation type, which may include small non-single-particle contributions that arise 
owing to the Ward-Takahashi identity. 

In QCD the relevant degrees of freedom change as the wavelength of the 
probe evolves. This feature is encoded, e.g., in the dressed-quark mass function, 
which is discussed in connection with Eq. ()A.18P . The nature of the mass function 
is model-independent and one consequence is that to a long wavelength probe a 
light-quark appears to have a large inertial mass ~ 350 MeV. 

Figure [2] expresses a nucleon current in which the primary degrees of freedom 
are dressed-quarks. Along with the Faddeev equation described in App. A, it is an 
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extension to baryons of the systematic and symmetry preserving rainbow-ladder 
truncation of QCD's DSEs, which provides a sound description of pseudoscalar 
and vector mesons and, in particular, a veracious description of the pion as 
both a Goldstone mode and a bound state of dressed-quarks [13]. It is a valid 
impulse approximation, which provides a systematically improvable continuum 
prediction for nucleon form factors. 

Subject to this understanding the question of whether a connection ex- 
ists between the spatial distribution of charge or magnetisation and the three- 
dimensional Fourier transform of a Sachs form factor involves a consideration of 
recoil-corrections experienced by dressed-quarks. The interpretation is appropri- 
ate if recoil corrections are small and can be calculated perturbatively. In that 
case the relevant expectation values in quantum mechanics are validly approxi- 
mated by the Fourier transform of the Sachs form factor. 

Consider the Breit frame and a photon probe with momentum Q = (0, 0, q, 0). 
In the scattering process this momentum is absorbed by the dressed-quarks 
within the proton. It is elastic scattering so all the dressed-quarks must recoil 
together, which means they can each be considered as absorbing a momentum 
fractioitl Q/2,. The magnitude of a recoil correction is then measured by the 
mass-squared scale 



Sr.- ^- 

We will consider that recoil corrections are small so long as 

1 



(11) 



Sr<^M\sr), (12) 

where M{s) is the dressed-quark mass function. This constraint means 

q< M(^) ^ g<0.4GeV, (13) 

a value determined from Eqs. (jA.lSp - ()A.22p . This momentum bound corre- 
sponds to a length-scale 

Ao.i = 0.49 fm = 0.57 rp, (14) 

where is the proton's charge radius. Hence in the three-dimensional Fourier 
transform of a Sachs form factor, recoil corrections are on the order of 10% or less 
throughout the domain r > 0.57 r^; namely, over 81% of the nucleons' volume. 
In measuring the total charge one must evaluate 

j-OO 

Qat = lim 47r / drr'^pN{r). (15) 

Ja 

It is interesting to reckon the amount of charge that is contained within the 
domain on which recoil corrections are not negligible. It is 

QV=4vr/ drr^pN{r). (16) 







■^Faddeev and Bethe-Salpeter amplitudes are peaked at zero relative momentum. Hence, the 
domain of greatest support in the impulse approximation calculation is that with each quark 



absorbing Q/2,. This is demonstrated explicitly, e.g., in Ref. [TT] . 
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A Gaufiian charge form factor can be used to obtain an algebraic and hence 
easily understood estimate; viz., 



G,(g2)=e-a^^^., (17) 



yields 



from which follows 



P.(0 = ^e--V[-|], (18) 



Q"/ = 0.19. (19) 

It is apparent that this region contains only 19% of the proton's charge. Expressed 
another way, the domain on which recoil corrections can be neglected contains 
81% of the proton's charge. (For the neutron's charge form factor the illustration 
can be made using a difference of two Gaufiians, each of which may be said to 
represent either the u- or d-quark contribution to the form factor.) If instead of 
Eq. (jl2p one were to consider recoil corrections as small for Sr < M'^{sr)/G, then 
the upper bound in Eq. ()16p is 0.48rp and the region contains only 12% of the 
proton's charge. 

On the other hand, recoil corrections are certainly large and essentially non- 
perturbative for 

Sr>M\sr) ^ g>lGeV, (20) 
a momentum boundary which corresponds to lengths 

Ai.o < 0.2 fm = 0.23 rp. (21) 

On this domain no quantum mechanical connection can be made between three- 
dimensional Fourier transforms of Sachs form factors and the density distribution 
of distinct charge and current carriers. It corresponds to 1.2% of the nucleon's 
volume and contains just 1.6% of the proton's charge. 

This analysis elucidates the circumstances under which the three-dimensional 
Fourier transform of a Breit-frame Sachs form factor can be viewed as providing 
a useful, qualitatively and semi-quantitatively reliable description of the config- 
uration space distribution of a nucleon's charge or magnetisation over dressed- 
quarks. Dressed-quarks are an emergent feature of QCD. The requisite conditions 
pertain within 81% - 99% of a nucleon's volume. Moreover, notwithstanding any 
caveats, Poincare invariant form factors are always a gauge of a hadron's struc- 
ture because they are a measurable and physical manifestation of the nature of 
the hadron's constituents and the dynamics that binds them together. 



5 Calculated Form Factors 



In the following two sections we present and discuss the results that our model 
of the dressed-quark core produces for nucleon form factors. Notably, we made 
significant modifications to the computer codes used to obtain the results in 
Ref. [31]. In addition to that described in App. D, which defines a convergent 
continuation of the Faddeev amplitude into the Breit frame, we succeeded in 
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Table 2. Probabilities associated with the F[ form factors evaluated at = 0, defined in 
App. E. 



0.474 


pP,c 

0.346 


pP,e 

0.180 


pP,s 

0.602 


pP,a 

0.254 


pp,m 

0.144 


0.441 


pP,c,u 

0.371 


0.188 


pp,s,li 

0.561 


0.294 


0.145 


pP,<l4 
0.345 


pp,c,d 

0.444 


pp,e,d 
0.210 


pp,s,d 

0.437 


pp.a,d 
0.414 


pp,m,d 

0.149 



reducing execution times by an order of magnitude. These two improvements 
enabled us to use a desk-top computer and obtain, within hours, numerically 
accurate results for the form factors on the domain G [0, 12] GeV^. 

In order to explain our results we must introduce our notation. The Pauli, 
Dirac and Sachs form factors are all represented by their usual symbols. Hence, 
the notation can be introduced by a single example. We choose the proton's Dirac 
form factor, Ff, and list the definitions in Appendix E: Form Factor Notation. 



It is also worth noting here that our analysis assumes rriu = m^. Hence 
the only difference between the u- and d-quarks is their electric charge. Our 
equations, computer codes and results therefore exhibit the following charge 
symmetry relations: 

CdFr = e^F^''' , e^Ff^" = e^F^ , i = 1,2. (22) 



6 Proton Form Factors 

6.1 Dirac proton 

In Fig. [3] we depict the proton's Dirac form factor and a breakdown into contri- 
butions from various subclasses of diagrams. The figures deserve careful study. 

The upper left panel shows the Q^-evolution of the quark, diquark and ex- 
change (or two-body) contributions to the form factor. Their = values mea- 
sure, respectively, the probability that the photon interacts with a bystander 
quark or a diquark correlation, or acts in association with diquark breakup: 

quark-Pf'^ = 0.47 : diquark-Pf '"^ = 0.35 : exchange- Pf''' = 0.18. (23) 

These and analogous probabilities are collected in Table [2j For Ff the diquark 
and exchange contributions switch in importance at ~ 3GeV^. Moreover, 
while the net result is always positive, the diquark contribution becomes negative 
at ~ 9GeV^. This panel here, and in kindred figures to follow, also displays 
a parametrisation of experimental results [H] for illustrative comparison with 
our computation. The manner by which that comparison should be understood 
is canvassed in Sect.O 

A radius can be associated with each of the form factors. We exemplify its 
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Figure 3. Proton's Dirac form factor. Upper left - Full result and decomposition according to 
diagram classes; lower left - flavour breakdown of these contributions, expressed in units of the 
magnitude of the relevant quark's charge; viz., |eu| — | or |ed| — |. Upper right - Full result 
and decomposition according to diagram diquark content; lower right - flavour breakdown of 
these contributions. A parametrisation of experimental data [44] is also presented in the upper 
left panel. A full explanation of the notation is provided in App. E. 



definition via i^f'"^; viz., 
and remark that 



d 



Q2=0 



(r?)2 = PI 



+pf'>n'+^{ 



1 



,p,e/ p,e\2 



1 



(24) 



(25) 



The calculated Dirac radii are reported in Table[3l Their values emphasise that so 
far as the Dirac form factor is concerned, the diquark component of the nucleon 
is softest. 

The lower left panel provides a flavour decomposition of the quark, diquark 
and exchange contributions to the form factor. While the other two u-quark 
components are positive definite, Ff'^'^ changes sign at Q2 ^ 9 GeV^. Up quarks 
are doubly represented in the proton and from Table [2] it is evident that they are 
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almost equally likely to be struck by a photon whether a bystander or a diquark 
participant. This explains the near equality of the radii associated with each term 
in the subclass of these form factor contributions in which a u-quark is struck. 
The same is not true for the d-quark, for which the probabilities show that it 
is more likely to be struck while a diquark participant. This signals that the 
d-quark is less free to move throughout the proton's volume and hence explains 
the small value of r^''^''^. 

The upper right panel of Fig. [3] shows the Q^-evolution of the contributions 
to Ff that involve a scalar diquark, an axial-vector diquark, or one of each. It is 
clear from Table [2] that the scalar diquark component of the proton is dominant. 
All contributions are positive definite, and the relative strength of the axial- 
vector and mixed contributions switches at ~ 5 GeV^ . From Table [3] one 
reads that the softest contribution to the proton's Dirac form factor is provided 
by diagrams involving an axial-vector diquark. One can picture this as stemming 
from the axial-vector correlation being more massive than the scalar and hence 
a bystander quark of any flavour ranges further from a collective centre-of-mass. 

The lower right panel provides a flavour decomposition of the diquark con- 
tributions just discussed. All u-quark components are positive definite. For the 
singly-represented d-quark, however, each of the form factors changes sign: F^''^'"^ 
becomes positive at ^ SGeV^; Ff'"''^ at _ SGeV^; and Ff'"''^ at 

~ 3 GeV^ . Axial- vector contributions to the Dirac form factor are the softest 
in each flavour sector. 

Evident in Table [3] is a notable feature of our calculation; viz., 

rf ^ > rf (26) 

Owing to charge symmetry this entails 

rf " > rf ^ (27) 

a result also obtained and explained in Ref. [23]. Equation ()26p follows from 
the presence of axial-vector diquark correlations in the nucleon. One reads from 
Table [2] that the proton's singly represented d-quark is more likely to be struck 
in association with an axial-vector diquark correlation than with a scalar, and 
form factor contributions involving an axial- vector diquark are soft. On the other 
hand, the doubly-represented u-quark is predominantly linked with harder scalar- 
diquark contributions. 



Table 3. Radii associated with F^, defined by analogy with Eq. (|24|l . All entries in fm. 
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rf 
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p.m 
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'1 
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' 1 
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rf 
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'1 

0.614 


p,s,u 

rf 
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rf ' 
0.656 


rf 

0.624 


'1 
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p.c.d 
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Figure 4. Proton's Pauli form factor. Left panel - Full result and decomposition according 
to diagram classes; right panel ~ full result and decomposition according to diagram diquark 
content. Form factors are expressed in magnetons defined by the calculated nucleon mass, Mn 
in Table [1] A parametrisation of experimental data [44] is also presented in the left panel. A 
full explanation of the notation is provided in App. E. 

6.2 Pauli proton 

In Figs. m and [5] we depict the proton's Pauli form factor and a breakdown into 
contributions from various subclasses of diagrams. 

The left panel of Fig.U] shows the Q^-evolution of the quark, diquark and 
exchange contributions to the form factor. Listed in Table HJ their = val- 
ues measure, respectively, the contribution to the proton's anomalous magnetic 
moment from the photon interacting with a bystander quark, a diquark or in 
association with diquark breakup. The net contribution from Diagrams 2 and 4 
in Fig. [2] is negative. This remains the case until ~ 9GeV^, at which point 
the net diquark contribution changes sign, as was also the case in Ff. The Pauli 
radii are listed in Table O from which it is evident that Diagrams 3, 5 and 6 in 
Fig. [2] provide the softest contribution. 



Table 4. Flavour and diagram breakdown of contributions to the proton's anomalous magnetic 
moment; viz., the form factors evaluated at = 0, measured in magnetons defined by the 
calculated nucleon mass, Mjv. 
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Figure 5. Proton's Pauli form factor: left panels - flavour breakdown of left panel in Fig.lH 
and right panels - flavour breakdown of right panel in Fig.|3] A full explanation of the notation 
is provided in App. E. 



The left panels in Fig. [5] provide a flavour decomposition of the quark, diquark 
and exchange contributions to the proton's Pauli form factor. We remark that 
is positive definite whereas changes sign at > 10 GeV^ and 

^ 17GeV^. (The latter should be interpreted qualitatively because our 
calculations are not truly reliable beyond 12 GeV^.) It is evident upon comparison 
between Tables [3] and [5] that the pattern exhibited by the Pauli radii is kindred 
to that of the Dirac radii, with the origin alike. 

The right panel of Fig.S] shows the Q^-evolution of the contributions to 
that involve a scalar diquark, an axial-vector diquark, or one of each. It is appar- 
ent from the figure and Table S] that diagrams involving the scalar correlation 
are dominant on a material domain. These contributions to a nucleon's Fad- 
deev amplitude have the simplest rest-frame spin-angular-momentum structure 
[21 [45]. We find that the scalar and axial- vector contributions are positive definite 
whereas the mixed contribution changes sign at Q2 > 8 GeV^. The latter provides 
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a larger contribution to the proton's magnetic moment than the axial-vector dia- 
gram. One reads from Table [5] that the softest contribution to the proton's Pauli 
form factor is provided by the axial-vector diquark diagrams. This was also the 
case for the Dirac form factor. However, in contrast to F^, the mixed contribu- 
tion to F2 is hardest, a result which owes primarily to Diagram 4 and the simple 
Ansatz we have made for the interaction therein; viz., Eq. ()C.35p . 

The right panels of Fig. [5] provide a flavour decomposition of the diquark 
contributions just discussed. It is curious that " < 0, a feature which highlights 
the presence and role of correlations in the nucleon's Faddeev amplitude. The 
associated form factor becomes positive at ~ 1.5 GeV^. The contribution 
with the simplest structure, Fg'*'", is positive definite whereas becomes 
negative at > lOGeV^. In association with the proton's d-quark, the axial- 
vector diagrams make a positive definite contribution, the scalar diquark form 
factor becomes negative at « 12 GeV^ and the mixed contribution is negative 
definite but small. 

It is apparent from Table [5] that 

r'/ > (28) 

which entails r2 > "^2 These orderings are the same as those exhibited by the 
Dirac radii, Eq. (I26p . but the separation in magnitudes is larger. The presence 
of axial-vector diquark correlations again plays a large role in producing these 
results. We note in addition that r2'" < r^'^ and r2''^ < r^''^, with the greater 
reduction in r2'". Indeed, it is almost uniformly true that the quark-core Pauli 
form factors are harder than their Dirac counterparts. The reduction is marked 
for 7-2'"'" and the only exception to the rule is rg''^''^. 

6.3 Pauli-Dirac proton ratio 

In Fig. [6] we plot a weighted ratio of Pauli to Dirac form factors; viz., 

""''^"^ ^ - (ln[QVi^])^ Fim ' " ' ^ - M^- ^'^^ 

A perturbative analysis that considers effects arising from both the proton's 
leading- and subleading- twist light-cone wave functions, the latter of which rep- 
resents quarks with one unit of orbital angular momentum, suggests that this 



Table 5. Radii associated with F^, defined by analogy with Eq. (|24|l . All entries in fm. 
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Figure 6. Solid curve - Calculated dimensionless ratio in Eq. (|29p with A = 0.44 and Mm 
in Table [1] Short-dashed line - constant at 0.21. Long dashed curve - Ratio evaluated with 
the experimental nucleon mass using the parametrisations in Ref. [33] . Boxes - Ratio evaluated 
with data from Ref. [47]; and Circles - from Ref. [48]. Dash-dot line - constant at 0.15. 



ratio should be constant for S> A'^, where yl is a mass-scale that is said to 
correspond to an upper-bound on the domain of soft momenta [l6] . 

We analysed our calculated result in this context and found that with A = 
0.44 this weighted ratio is a constant = 0.21 on > 4.3; by which we mean 
that the rms relative error with respect to the straight-line fit is 0.34% with an 
associated standard deviation of 0.24%. These numbers increase as the minimum 
value of included in the fit is decreased and, moreover, the value of A comes 
to depend on this minimum value. 

In the figure we also plot the ratio in Eq. (j29p as evaluated from extant exper- 
imental data, available on the domain E [3.9, 6.3]. Using A = 0.44, the result 
is described by a constant = 0.15 with rms relative error 1.5% and an associated 
standard deviation of 0.98%. It is evident in the figure that on the domain for 
which the ratio is well described by a constant, our model produces a result that 
lies above the experimental data. This is because thereupon our calculation un- 
derestimates experimental results for Ff by ~ 15% and overestimates those for 
F2 by a similar amount. (See Sect. [8] for details.) 

It is curious that what might appear to be a low mass-scale, A = 0.44M7V, 
should serve to produce a constant value for the ratio in Eq. (j29p |31j . In seeking 
to understand the origin of this scale we analysed the pointwise behaviour of our 
calculated Faddeev amplitude. The dominant functions are si, ^3, A^, which was 
to be expected given the associated Dirac structures (see Eqs. (jA.SP - (jA.lOp .) 
A Gaufiian can be fitted to the leading Chebyshev moment of each of these 
functions. That procedure yields the following widths (in units of M^r): 



= 0.48 , w^i = 0.47 , uj^i = 0.46 . 



(30) 
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Table 6. Radii associated with G^, defined by analogy with Eq. (|24|l . NB. The value in this 
table yields Mjvr^ = 4.01 cf. experiment [49] Mjvr|; = 4.18. Tabulated entries in fm. 
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The similarity between these widths and A is notable. It highlights the point 
that while A per se is not an elemental input to our calculation, such a mass- 
scale can arise dynamically as a derivative quantity which may be expressed in 
the relative-momentum support of the Faddeev amplitude. A challenge now is 
to determine whether an algebraic relationship exists between A in Eq. (|29p and 
the widths characterising the Faddeev amplitude. 

6.4 Sachs proton electric 

In Fig. [7] we present the proton's Sachs electric form factor and a separation 
into contributions from various subclasses of diagrams. While in principle, given 
Eqs. (HI), these panels contain no information that cannot be constructed from 
material already presented, they are nevertheless practically useful and informa- 
tive. 

The upper left panel shows the Q^-evolution of the quark, diquark and ex- 
change (or two-body) contributions to G^((5^). Their = values have pre- 
cisely the same value and interpretation as those associated with the Dirac form 
factor, which are presented in Table [2l It is notable that the quark contribu- 
tion; namely. Diagram 1 in Fig. [21 possesses a zero at ~ 3.0 GeV^. It is only 
because the diquark contribution remains positive until ~ 9.0 GeV^ and the 
exchange contribution is positive definite that the complete result for G^{Q'^) 
does not exhibit a zero until Q2 ~ s.OGeV^I We list the Sachs radii in Table [H 
In comparison with the Dirac radii in Table [3l they are relatively uniform owing 
to Foldy-term contributions. 

The lower left panel provides a flavour decomposition of the quark, diquark 
and exchange contributions to the form factor. Once more their = values 
have precisely the same value and interpretation as those associated with the 
Dirac form factor. G^" has a zero at ~ 9.0 GeV^ and no contribution to 
G^" is positive definite: G^^'" possesses a zero at ^ 3.0 GeV^; G^^'^'^ at 
Q2 ~ 10.0 GeV^; and G^''''' is negative on the domain 4.0 < < 7.0GeV2. On 
the other hand, G^'^ has a zero at ~ 10.0 GeV^ but G^''''^ is negative definite. 

zero in G^{Q^) was seen in the light-front constituent-quark model of Ref. [50) . In Ref. [51] it 
was shown to be a property of the scalar-diquark Faddeev model of Ref. [37] but its appearance 
and location were argued to be dependent on dynamics, consistent with Refs. [521 153) and the 
present study. 
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Figure 7. Proton's Sachs electric form factor. Upper left - Full result and decomposition ac- 
cording to diagram classes; lower left - flavour breakdown of these contributions, expressed 
in units of the magnitude of the relevant quark's charge. Upper right - Full result and de- 
composition according to diagram diquark content; lower right - flavour breakdown of these 
contributions. A parametrisation of experimental data 31] is also presented in the upper left 
panel. A full explanation of the notation is provided in App. E. 



G^"''^ possesses a zero at « 11.0 GeV^ and G^"''^ at « 2.0 GeV^. We list 
the u- and d-quark Sachs radii in Table [6l Their values are readily computed and 
understood from Tables [3] and HI 

The upper right panel of Fig.[7]shows the Q^-evolution of the contributions to 
that involve a scalar diquark, an axial-vector diquark, or one of each. has 
a zero at « 10.0 GeV^ and G^/ at w 3.0 GeV^, whereas G^"" is positive 
definite. The lower right panel provides a fiavour decomposition of the diquark 
contributions. G^"'" exhibits a zero at ^ 10.0 GeV^, G^^"'" at ~ 5.0 GeV^ 
and G^ ' is negative on the domain 2.0 < < e.OGeV^. On the other hand, 
G^/''^ passes through zero at « 11.0 GeV^ and G^""'"^ at w 2.0 GeV^ but 
Gg*^'*^ is negative definite. The associated Sachs radii are listed in Table [6l 
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Figure 8. Proton's Sachs magnetic form factor. Left panel - Full result and decomposition 
according to diagram classes; and right panel - full result and decomposition according to 
diagram diquark content. Form factors are expressed in magnetons defined by the calculated 
nucleon mass, Mjv in Table [T] A parametrisation of experimental data [U is also presented in 
the upper left panel. A full explanation of the notation is provided in App. E. 

6.5 Sachs proton magnetic 

In Figs. [8] and [9] we depict the proton's Sachs magnetic form factor and a sep- 
aration into contributions from various subclasses of diagrams. Again, while in 
principle these panels only contain information that can be constructed from ma- 
terial already presented, they are nonetheless practically useful and informative. 

The left panel of Fig. [8] shows the Q^'^volution of the quark, diquark and 
exchange contributions to the form factor. G^, and are positive definite 
and monotonically decreasing. On the other hand, the net contribution from 
Diagrams 2 and 4 in Fig. [21 namely, G^f, is uniformly small, negative in the 
vicinity of ^ 0.5 GeV^ and again for > 8 GeV^. The pattern is qualitatively 
similar in the flavour breakdown of these form factors, depicted in the left panels 
of Fig.El 



Table 7. Flavour and diagram decomposition of contributions to the proton's magnetic mo- 
ment; viz., the form factors evaluated at = 0, measured in magnetons defined by the 
calculated nucleon mass, Mjv. Experimentally [49], A^p — 2.79. 
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Figure 9. Proton's Sachs magnetic form factor. Left panels - Flavour breakdown of contribu- 
tions in left panel of Fig.|8l Right panels - Flavour breakdown of contributions in right panel of 
Fig-El A full explanation of the notation is provided in App. E. 



The right panel of Fig. [8] exhibits the Q^-evolution of the contributions to 
G^j^j that involve a scalar diquark, an axial-vector diquark, or one of each. All 
contributions are positive definite, diagrams involving only a scalar diquark are 
dominant and contributions involving at least one axial-vector diquark are uni- 
formly of comparable magnitude. The flavour breakdown is contained in the 
right panels of Fig.[9l all contributions are positive definite except G^^'"^, which 
is uniformly small but becomes positive at ~ 9.0 GeV^ and remains so un- 
til 17GeV^. (NB. The latter should be interpreted qualitatively because 
the feature appears at a larger value of than we consider our computation 
reliable.) 

In Table [7] we list the = values of all the form factors, which measure, 
respectively, the contributions to the proton's magnetic moment. These values 
can be obtained from = Ff(0)P" -|- Kp, where, e.g., a = {p,q),{p,c), etc. 
The magnetic radii are listed in Table [8l The general pattern has electric radii 
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Figure 10. Result for the normalised ratio of Sachs electric and magnetic form factors computed 
with four different diquark radii, rj^+. Data: diamonds - [3H]i squares - triangles - [56j : 
and circles [57l. 



exceeding magnetic radii. The few exceptions are easily explained. For example, 
r^"'" < r^'^: this is primarily because F2'"''^{Q'^) < and of significant magni- 
tude in the neighbourhood of = 0. As already noted, it is curious that this 
contribution to the proton's anomalous magnetic moment is negative. 

6.6 Sachs electric-magnetic proton ratio 

We plot fip G^^iQ"^) / G\^{Q'^) in Fig.[TO]in comparison with contemporary data. A 
sensitivity to the proton's electromagnetic current is evident, here expressed via 
the diquarks' radius. Irrespective of that radius, however, the proton's electric 
form factor possesses a zero and the magnetic form factor is positive definite. 
On < 3 GeV^ our result lies below experiment. As discussed in Sect. [HI this 
can likely be attributed to our omission of so-called pseudoscalar-meson-cloud 



Table 8. Radii associated with G\j, defined by analogy with Eq. (|24|) . NB. The value in this 
table yields Mnr^^j ~ 3.23 cf. experiment [54] Mnv'^j = 3.99. Tabulated entries in fm. 
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Figure 11. Neutron's Dirac form factor. Left panel - Full result and decomposition according 
to diagram classes; Right panel - Full result and decomposition according to diagram diquark 
content. A parametrisation of experimental data [44] is also presented in the left panel. A full 
explanation of the notation is provided in App. E. 



contributions. 



7 Neutron Form Factors 

7.1 Dirac neutron 

In Fig.[TT] we depict the neutron's Dirac form factor and a decomposition into 
contributions from various subclasses of diagrams. Owing to charge symmetry, 
Eqs. ()22p . it is unnecessary to present a flavour breakdown. For example, with the 
normalisation used in our figures, the curve that would be denoted by F^'^{Q'^) 
is simply negative-i^f''^((5^) drawn from Fig.[3l 

In addition to that of FJ^ itself, the left panel panel depicts the Q^-evolution 
of the quark, diquark and exchange contributions to this form factor. F"' and F^''^ 
are negative definite, and i^"'^ is only positive for < 0.5 GeV^. On the other 
hand, the diquark contribution; viz., F^'"^, is positive until ~ 12GeV^. The 
right panel renders the Q^-dependence of contributions from diagrams containing 
a scalar diquark, an axial-vector diquark or one of each. F^'^ is negative definite 
and is negative for > 2 GeV^. F"'"" is smah at = (only 3% of the 
other two form factors) and negative for > 0.1 GeV^. These features reflect: 
the dominant role played in the Faddeev amplitude by the positively-charged 
[ud] scalar diquark; the fact that the u-quark is singly-represented and only a 
bystander in combination with an axial-vector diquark; and the softness of the 
diquark correlations, which ensures that only a bystander quark can participate 
in the scattering process at large-Q^. 

We list computed Dirac radii connected with the neutron in Table [H Two 
entries are imaginary because the associated form factors have an inflexion point 
away from = 0. We do not currently attribute any real signiflcance to this 
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Figure 12. Neutron's Pauli form factor. Left panel - Full result and decomposition according 
to diagram classes; Right panel - Full result and decomposition according to diagram diquark 
content. A parametrisation of experimental data [44] is also presented in the left panel. A full 
explanation of the notation is provided in App. E. 

local feature, which for the neutron is particularly sensitive to details of the 
Ansatz employed for Diagrams 5 and 6 in Fig.[2l namely, the two-body piece of 
the current which is not yet well constrained. 

7.2 Pauli neutron 

In Fig. [12] we depict the neutron's Pauli form factor and a decomposition into 
contributions from various subclasses of diagrams. Once more, owing to charge 
symmetry, Eqs. (j22p . it is unnecessary to present a flavour breakdown. For exam- 
ple, with the normalisation used in our figures, the curve that would be denoted 
by is simply negative-F2'^((5^) drawn from Fig. [J] 

The left panel panel depicts the Q^-evolution of F2 itself, and that of the 
quark, diquark and exchange contributions to this form factor. F2, F^'*^ and F^'*^ 
are negative definite on the domain within which we consider our calculations 
accurate, and F^''^ is negative until ~ 12GeV^. The right panel portrays the 
Q^-dependence of contributions from diagrams containing a scalar diquark, an 
axial-vector diquark or one of each. F^'* and F^'" are negative definite, and ^2'"^ 

Table 9. Radii associated with Fi , defined by analogy with Eq. (|24p except when the form 
factor vanishes at — 0, in which case = —6F'(Q^ — 0). An imaginary result signifies a 
negative rms radius. This convention enables a straightforward comparison between the length- 
scale associated with different radii. All entries in fm. 
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y=Q^/M^ 

Figure 13. Solid circles and solid curve - Dimensionless ratio in Eq. (|29|) calculated for the 
neutron, with A = 0.44 and Afjv in Table [1] Dashed curve - Right-hand-side of Eq. (|3ip . 
Experimental results: down triangles - Ref. [58] . 



is negative for < 5 GeV^ and always small in magnitude. These features are 
consistent with those of the Dirac form factor. 

We list computed anomalous magnetic moments and Pauli radii connected 
with the neutron in Table [101 The small value of may be understood via a 
cancellation between d{ud)i+ and u{dd)i+ contributions. Along with the small 
value of kt, Eq. (|C.36h . this explains the size of k™. With the exception of the 
uniformly small -F^''^, the Pauli radii follow the same pattern as those of the 
proton. 

7.3 Neutron Pauli-Dirac neutron ratio 

In Fig. [13] we plot the weighted ratio of Pauli to Dirac form factors in Eq. (129p 
for the neutron. This ratio is constant for the proton, Fig.[6l however, that is 
not the case for the neutron. Moreover, with our calculated neutron form factors 



Table 10. Upper rows: Diagram decomposition of contributions to the neutron's anomalous 
magnetic moment; viz., the form factors evaluated at = 0, measured in magnetons 
defined by the calculated nucleon mass, Afjv. Experimentally [49], /in = —1.91. Lower rows: 
Radii associated with F^, defined by analogy with Eq. (|24|l . These entries in fm. An imaginary 
result signifies a negative rms radius. 
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Figure 14. Neutron's Sachs electric form factor. Left panel - Full result and decomposition 
according to diagram classes. Right panel - Full result and decomposition according to diagram 
diquark content. A parametrisation of experimental data [44] is also presented in the left panel. 
A complete explanation of the notation is provided in App. E. 



there is no value of A for which this ratio assumes a constant value. 

The apparent cause of this behaviour is a zero in F2{Q'^) at ~ 18GeV^. 
This point lies beyond the upper bound of the domain within which we con- 
sider our computation reliable. On the other hand, its presence does influence 
the evolution of the ratio. This can be seen by analysing the ratio using Fade 
approximants on subdomains of G [4, 12]GeV^, which consistently yields a 
best fit that possesses a zero at « 18GeV^; e.g., 

Rn (n2^ F^iQ^) ^ 2.85 + 0.274 _ 0.0409 q4 

^ (lnQ2/yi2)2 ^n(g2) _1 + 1.93Q2 

It seems therefore that the zero is not simply the result of inaccurate numerical 
analysis but either a property of the model itself or an artefact of the numerical 
method; namely, the Chebyshev expansion described in Appendix D: Chebyshev 
expansion. We are working to resolve this issue. 



7.4 Sachs neutron electric 

In Fig. [13] we present the neutron's Sachs electric form factor and a separation 
into contributions from various subclasses of diagrams. Once more, owing to 
charge symmetry, Eqs. (j22|) . it is unnecessary to present a flavour breakdown. 
For example, with the normalisation used in our figures, the curve that would 
be denoted by G^'^{Q'^) is simply negative-G^"(Q2) drawn from Fig.[71 

In addition to that of itself, the left panel panel depicts the Q^.gyohition of 
the quark, diquark and exchange contributions to this form factor. Each exhibits 
a zero, with that for the net result lying at ~ llGeV^. In the right panel 
we plot the Q^-dependence of contributions from diagrams containing a scalar 
diquark, an axial-vector diquark or one of each. G^'^ is positive on the domain 



I. C. Cloet, et al. 



25 



g2 G [0.1,ll]GeV2 and G^'™ is negative for > iGeV^, whereas G^'" is 
positive definite. 

These features are again consistent with intuition. For example, the behaviour 
of G^'^ . It is negative at smah-Q^ because the scalar diquark component of the 
Faddeev amplitude is dominant and that is paired with a d-quark bystander in 
the neutron. This dressed-quark is responsible for the preponderance of negative 
charge at long range. G^'^ is positive at large because F2 dominates on that 
domain, which focuses attention on the axial-vector diquark component of the 
Faddeev amplitude. The positively charged u-quark is most likely the bystander 
quark in these circumstances. 

Another interesting illustrative case is provided by G^'", which is positive 
definite. As already noted, the n-quark is the most probable bystander in the 
neutron's axial-vector diquark configuration and this explains the preponderance 
of positive charge at small Q^. This plus the fact that the current's only hard 
component is that involving a bystander quark also explains the positive charge 
at large Q^. The form factor remains positive in the intermediate region because 
the term which could interfere; viz., d[ud]Q+, involves individual charges with 
smaller magnitude. 

We list computed Dirac radii connected with the neutron in Table [H Two 
entries are imaginary because the associated form factors have an inflexion point 
away from = 0. As just explained, the origin of such behaviour lies in in- 
terference, mediated by the current, between components in the incoming and 
outgoing neutrons' Faddeev amplitudes. 

7.5 Sachs neutron magnetic 

In Fig. [15] we present the neutron's Sachs magnetic form factor and a decompo- 
sition into contributions from various subclasses of diagrams. Again, owing to 
charge symmetry, Eqs. (|22p . it is unnecessary to present a flavour breakdown. 
For example, with the normalisation used in our figures, the curve that would 
be denoted by G^{Q'^) is simply negative-G^j'^((5^) drawn from Fig.O 

In the left panel panel we draw the Q^-evolution of G^ itself, and that of the 
quark, diquark and exchange contributions to this form factor. G'^j, G^^ and 

Table 11. Upper rows - Radii associated with G^, defined by analogy with Eq. (|24|l except 
when the form factor vanishes at — 0, in which case = -6F'(Q2 = 0). An imaginary 
result signifies a negative rms radius. NB. The value in this table yields M^{r^)^ — —(1.36)^ 
cf. experiment [49] M^{r^)^ = —(1.62)^. Lower rows - Radii associated with G^/, defined by 
analogy with Eq. ((24l): MNr^' = 3.17 cf. experiment [54] M„r^ = 4.24. Tabulated entries in 
fm. 
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Figure 15. Neutron's Sachs magnetic form factor. Left panel - Full result and decomposition 
according to diagram classes. Right panel - Full result and decomposition according to diagram 
diquark content. A parametrisation of experimental data [4j is also presented in the left panel. 
A complete explanation of the notation is provided in App. E. 

are negative definite. On the other hand, is uniformly small, owing to 
cancellations between and Fl*. It begins negative, is positive in the vicinity 
of = 0.5 GeV^ and again for > lOGeV^. The right panel portrays the 
Q^-dependence of contributions from diagrams containing a scalar diquark, an 
axial-vector diquark or one of each. All are negative definite. 

We list the computed magnetic radii connected with the neutron in Table [TTl 
The magnetic moments are the same as the anomalous moments in Table [TOl 
With the exception of G^^, which at small is roughly a factor of five smaller 
than G'f^^, the neutron radii follow the same pattern as those of the proton. 

7.6 Sachs electric-magnetic neutron ratio 

We plot fJ'nG^{Q'^) /G'^^{Q'^) in Fig. [TBI The figure illustrates a quantitative sen- 
sitivity of our results to the neutron's electromagnetic current, here expressed via 
the diquarks' radius. Notwithstanding this, the qualitative features are robust, 
with G^{Q'^) possessing a zero at > 10 GeV^. In contrast to the behaviour in 
Fig. llOl here the zero moves to smaller with increasing diquark radius. The 
effect of our omission of meson cloud contributions is again evident at small . 

8 Chiral Corrections 

The framework we have described hitherto defines a dressed-quark core contri- 
bution to the nucleons' electromagnetic form factors. As with the mass |281 129|. 
the nucleons' magnetic moments, and charge and magnetic radii receive mate- 
rial contributions from the so-called pseudoscalar meson cloud |59^ 160] . There 
are two types of contribution: regularisation-scheme-dependent terms, which are 
analytic functions of m in the neighbourhood of vanishing current-quark mass. 
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Figure 16. Result for the normalised ratio of Sachs electric and magnetic form factors for 
the neutron computed with two different diquark radii. Short-dashed curve: parametrisation of 
Ref. [44]. Down triangles: data from Ref. [58] . 



m = 0; and nonanalytic scheme- independent terms. 

For magnetic moments and radii the leading-order scheme-independent con- 
tributions are [M] 

KH-n/plNA — ^ ^^2 p '""^ ' 

{rn/p)NA - (33) 



^^^^^^^ " s^^'^Mf^ + ie^;;^:' ^^^^ 



where, experimentahy, qa = 1.26, /tt = 0.0924 GeV= 1/(2.13 fm) and fiy = 
fJ'P — /^n. = 4.7. These terms reduce the magnitude of both neutron and proton 
magnetic moments, and increase the magnitudes of the radii. 

Whilst these scheme-independent terms are important, at physical values of 
the pseudoscalar meson masses they do not usually provide the dominant contri- 
bution to observables. That arises from the regularisation-parameter-dependent 
terms, as is apparent for baryon masses in Ref. [28] and for the pion charge radius 
in Ref. [62]. This is particularly significant for the neutron's charge radius [31] 
and for the magnetic moments, in which connection the regularisation-scheme- 
dependent terms provide a nonzero contribution in the chiral limit and have the 
net effect of increasing \fiN\- 

Owing to the importance of the chiral loops' regularisation-parameter- 
dependent parts we estimate the corrections using modified formulae, which 
incorporate a single parameter that mimics the effect of regularising the inte- 
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Table 12. Quark-core and pseudoscalar meson loop [Eqs. (|35p - (|37p ] contributions to the mo- 
ments and radii, calculated at the physical current-quark mass, Eq. (|X23)) . The radi i are listed 
in fm'^. Experimental values are quoted from Ref. [4^, where available, and otherwise from 
Ref. \M- 
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wherein /i^^^ are the chiral hmit values of the meson loop contributions and 
A = 0.3GeV= l/[0.66fm] is a regularisation mass-scale. NB. As required phys- 
ically, the loop contributions vanish when the meson mass is much larger than 
the regularisation scale: very massive states must decouple from low-energy phe- 
nomena. 

In Table [T2] we exemplify the effect of the corrections in Eqs. (j35p - (|37p to nu- 
cleon static properties. The quark-core values are collected from Tables EHSl and 
[TOHl II herein. The sensitivity of the neutron's charge radius is apparent. In rela- 
tion to the magnetic moments, a recent estimate from numerical simulations of 
lattice-regularised QCD |64j gives the following chiral-loop contributions to the 
nucleons' magnetic moments at the physical pion mass: ji!^ = —0.40, ^u^ = 0.24, 
which are obtained with /i^^ = —1.05, = 0.88 in Eq. (j35p . These results 
in conjunction with the experimental values point to quark-core magnetic mo- 
ments of n^'^'^^ = —1.51, /ip^^'^'' = 2.55, which compare well with our computed 
moments □ 

It is plain in Table [T2] that pseudoscalar meson loops alter the proton's mag- 
netic radius more than its electric radius. Indeed, without fine tuning, these 
two initially rather different radii are brought into agreement. As observed in 



*The magnetic moment values in Row 2 of the Table differ slightly (< 8%) in magnitude from 
those reported in Ref. [33] because an extrapolation is necessary to obtain Gm{0) and herein 
we've used a [0, 2] Fade as opposed to a simple quadratic. Were this significant, it could be 
corrected by a minor (~ 10%) adjustment of fii+ and kt- 
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Figure 17. DifFerence between our calculated Pauli form factor and the parametrisation of 
experimental data in Ref. [44], each normalised by the appropriate anomalous magnetic moment 
at = 0: dashed curve ~ proton; solid curve - neutron. The for which the difference reaches 
20% of its peak value is indicated in each case by a vertical dotted line. 



Ref. [31], this is important in relation to Fig.[TU] because it explains why the 
quark core result disagrees with data at small momentum transfers. Namely, in 
the neighbourhood of = one has 

and so with > r^, as is the case for the quark core contribution, the ratio 
falls immediately with increasing Q^. This is the behaviour in Fig.llUi However, 
experimentally, and with addition of a pseudoscalar meson cloud to our quark 
core, rp = Vp. Therefore the complete ratio varies little on < < 0.6 GeV^. 

The analysis in this section is rudimentary. Nonetheless it illustrates that the 
dressed-quark core defined by our Faddeev equation is uniformly compatible with 
augmentation by a sensibly regulated pseudoscalar meson cloud. We emphasise 
that by construction our Faddeev equation explicitly excludes all diagrams that 
can be associated with that cloud and so a question of overcounting cannot arise. 

It is nevertheless reasonable to inquire into the domain of momentum trans- 
fer upon which pseudoscalar meson loops can contribute materially to nucleon 
form factors. Regarding this it is relevant to observe that in a meson- nucleon 
coupled-channel analysis of the A transition form factors the cloud con- 

tributes 50% of the Ml form factor's magnitude at = but is insignificant 
by ~ 2M^ [65]. We address this question via Fig. 1171 which compares our 



'TV 

computed dressed-quark core Pauli form factors with a contemporary parametri- 
sation of experimental data [44]. The differences depicted are consistent with 
loop corrections providing a necessary quantitative contribution that is impor- 
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tant until w 2-3A/^. An analogous figure for the Dirac form factors presents 
a comparable picture, although the differences are an order of magnitude smaller 
and have longer tails. 

9 Epilogue 

We described a calculation of a dressed-quark core contribution to nucleon elec- 
tromagnetic form factors. This core is defined by the solution of a Poincare covari- 
ant Faddeev equation in which dressed-quarks provide the elementary degree of 
freedom and quark-quark correlations are formed therefrom. The two parameters 
in the Faddeev equation are diquark masses. They are set by fitting to required 
nucleon and A masses. We allowed one parameter in the nucleon-photon ver- 
tex; viz., the diquark charge radius. Contemporary continuum calculations and 
comparison with extant data indicate that this radius should be commensurate 
with the pion's charge radius. From this foundation we provided a comprehensive 
analysis and explanation of the form factors. 

A feature of our study is the separation of form factor contributions into 
those from different diagram types and correlation sectors, and subsequently a 
flavour separation for each of these. In this way we obtained, for example, Tabled 
which shows amongst other things that the probability of the photon striking a 
bystander quark in the proton is 47%. It also enables us to determine, Eq. (f27l) . 
that r"'" > r"''^; i.e., that the neutron's n-quark Dirac radius is greater than that 
of the d-quark, and explain the result in terms of the presence of axial-vector 
diquark correlations. The dressed-quark magnetic radii have the same ordering. 

From our extensive body of results we will here highlight just a few more. 
For the proton a weighted ratio of Pauli to Dirac form factors is constant on a 
domain that begins at Q^/M^ ~ 4, Fig.[6l We correlated this behaviour with 
the momentum space width of the dominant elements in the proton's Faddeev 
amplitude, Eq. (j30p . On the other hand, the same ratio for the neutron is not 
constant on any domain accessible in our calculation. Fig. [131 In addition, the 
ratio of Sachs electric and magnetic form factors for the proton exhibits a zero, 
Fig.Hni Its position depends on correlations in the Faddeev amplitude and details 
of the nucleon-photon current. Our current best estimate for the zero's location 
is ~ 8GeV^. A similar ratio for the neutron passes through zero at ~ 
llGeV^, Fig.dSl 

We have defined the nucleon's dressed-quark core via a Poincare covariant 
Faddeev equation and have seen that pseudoscalar meson loops can be added 
in a sensible fashion. The framework is successful and instructive, and unifies a 
phenomenological description of mesons with that of baryons. Yet it is simple 
enough to allow access to numerous form factors and large values of momentum 
transfer. Importantly, our approach enables one to chart the interplay between 
the firmly established and material momentum-dependent dressing of QCD's 
elementary excitations and the observable properties of hadrons. In the near 
term it will be applied to nucleon excited states and transition form factors 
so as to elucidate their dependence on these fundamental features of QCD. A 
medium term goal is to extend Ref. [23] and provide a simultaneous description 
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of meson and baryon observables using a single interaction in a truncation of 
QCD's Dyson-Schwinger equations that can systematically be improved. 
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Appendix A: Faddeev Equation 

A.l General structure 

The nucleon is represented by a Faddeev amplitude 

i^/ = l^/j + 0/2 + if's , (A.l) 

where the subscript identifies the bystander quark and, e.g., !f'i,2 are obtained from •Z'a by a 
cyclic permutation of all the quark labels. We employ the simplest realistic representation of 
i/. The spin- and isospin-1/2 nucleon is a sum of scalar and axial- vector diquark correlations: 

<p3{p^,ai,n)=jVt +JVt, (A.2) 

with {pi,ai,Ti) the momentum, spin and isospin labels of the quarks constituting the bound 
state, and P = pi + P2 + Ps the system's total momentum. 
The scalar diquark piece in Eq. (|A.2p is 

AAf (p„ a.,n) = [r°^ip[i2] ; K)]IY^^ {K) [S{t, P)u{P)]l\ , (A.3) 

where: the spinor satisfies ( [Appendix 'S>:\ Euclidean Conventions) 

(i-y- P + M) u{P) = = u{P) (ij- P + M) , (A.4) 

with M the mass obtained by solving the Faddeev equation, and it is also a spinor in isospin 
space with ip+ = col(l,0) for the proton and (fi- = col(0, 1) for the neutron; K — pi + P2 ~: 
P{i2}, P[i2] = Pi ^ P2, £ ■= (— P{i2} + 2p3)/3; Zi"^ is a pseudoparticle propagator for the scalar 
diquark formed from quarks 1 and 2, and _r°^ is a Bethe-Salpeter-like amplitude describing their 
relative momentum correlation; and 5, a 4x4 Dirac matrix, describes the relative quark-diquark 
momentum correlation. (5, F^' and zi" are discussed in Sect. IA.2l ') The colour antisymmetry 
of 'Fs is implicit in F , with the Levi-Civita tensor, ecic2C3i expressed via the antisymmetric 
Gell-Mann matrices; viz., defining 

{f^'=iA^//^ = -^A^//^ = iA'}, (A.5) 

then ecic2C3 = {H''^)c,c2- [See Eqs. ^EM, ^EM-] 
The axial-vector component in Eq. HA.2|I is 

M'\p,,a,,n) = [t'F'^''(^p[,2];K)]l\^^^AlUK)[Al{£-,P)u{P)]l\, (A.6) 

where the symmetric isospin-triplet matrices are 

t+ = + r') , t" ^r\t- = ±{r°~ r') , (A.7) 

and the other elements in Eq. (jA.6[l are straightforward generalisations of those in Eq. (|A.3|l . 

The general forms of the matrices S{£;P) and Al{£;P), which describe the momentum 
space correlation between the quark and diquark in the nucleon are described in Refs. [21 145] . 
The requirement that S{£;P) represent a positive energy nucleon entails 



S{£;P)^si{£;P)Ii,+ (ij-i-i-PlD^j S2{£;P), 



(A.8) 
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where (I-D)rs = <5rs, i = 1, P = —1. In the nucleon rest frame, si,2 describe, respectively, the 
upper, lower component of the bound-state nucleon's spinor. Placing the same constraint on 
the axial- vector component, one has 



(A.9) 



where {e^ ^ + £ ■ P P,, ■y^ ^ -f, + -/ ■ P P^) 



= 7 • £^ P, . 



At 



'iPu : 



One can now write the Faddeev equation satisfied by $'3 as 

d^£ 



-f£^£^ 



At 



S{k-P)u{P) 
AUk;P)n{P) 



The kernel in Eq. ()A.11|I is 



M{k,£- P) 



(2^) 

Moo 
(Mio)] 



■Mik,£; P) 



S{£;P)u{P) 
AU£; P)u{P) 



{Mc 



with 



Moo = r^^ikg - £gg/2- £gg) S' {£gg - kg) P" ' {£g - fc,,/2; -kgg) S{£g) A" ' {£gg) 



0+, 



(A.IO) 



(A.ll) 



(A.12) 



(A.13) 



where: £g=£ + P/3, kg = k + P/3, £gg 
denotes matrix transpose; and 



-^-^2P/3, kg 



-k + 2P/3 and the superscript "T" 



{M0l)i=t' r^'^{kg-£gg/2;£gg) 



{^qq kq) ^ i^q kqq/l] k qq) S {£ q) /\ {£ qq) ^ 



{Ml0)l^r°\kg-£gg/2-£gg) 

X S i^^qq ~ kq)t, (^£q — kqq/2] ^gg) S(^iq) /\ ('^99 ) 5 

iMll)%=e ryikg-£gg/2;£gg) 

X S {£qq — kq) t [£q — kqq/2] —kqq) S{iq) i^qq) • 



(A.14) 
(A.15) 
(A.16) 



A. 2 Kernel of the Faddeev equation 

To complete the Faddeev equations, Eq. HA.11|) . one must specify the dressed-quark propagator, 
the diquark Bethe-Salpeter amplitudes and the diquark propagators. 



A. 2.1 Dressed-quark propagator 

The dressed-quark propagator has the general form 



Sip) 



-i-y -pavij) ) + crs{p ) = l/[i7 ■ P^(p ) + B{p )] 



(A.17) 



and can be obtained from QCD's gap equation; namely, the Dyson-Schwinger equation for the 
dressed-quark self-energy [5] • It is a longstanding prediction of DSE studies in QCD that for 
light-quarks the wave function renormalisation and dressed-quark mass: 



Zip") = l/Aip") , M(p^) = B{p^)/A{p') , 



(A.18) 



respectively, receive strong momentum-dependent corrections at infrared momenta O |4j [5]: 
Z(p'^) is suppressed and M{p^) enhanced. These features are an expression of dynamical chiral 
symmetry breaking (DCSB) and, plausibly, of confinement [14]. The enhancement of M{p^) 
is central to the appearance of a constituent-quark mass-scale and an existential prerequisite 
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for Goldstone modes. These DSE predictions are confirmed in numerical simulations of lattice- 
regularised QCD [6], and the conditions have been explored under which pointwise agreement 
between DSE results and lattice simulations may be obtained [71 [8l 166). 

The impact of this infrared dressing on hadron phenomena has long been emphasised [11] 
and, while numerical solutions of the quark DSE are now readily obtained, the utility of an 
algebraic form for S{p) when calculations require the evaluation of numerous multidimensional 
integrals is self-evident. An efficacious parametrisation of S{p), which exhibits the features 
described above, has been used extensively in hadron studies [BT]- It is expressed via 

J^{2{x + m')) + J^ibix) T{b:ix) [bo + biTiex)] , (A.19) 

av{x) = _!_ [1 - jr(2(a; + ^2))] , (A.20) 
X + m-' 

with X = p^/A^, fh — m/X, 

T{x) = , (A.21) 

X 

o'six) — Xas{p'^) and av{x) — A^(Jv(p^). The mass-scale, A — 0.566 GeV, and parameter 
valuefl 

^ bo bi b2 bz 

0.00897 0.131 2.90 0.603 0.185 ' y ■ I 

were fixed in a least-squares fit to light-meson observables [681 169[ . The dimensionless u = d 
current-quark mass in Eq. (|A.22p corresponds to 

m = 5.08 MeV =: m^^^'' . (A.23) 

The parametrisation yields a Euclidean constituent-quark mass 

Mi^d = 0.33 GeV, (A.24) 

defined as the solution of — M^{p'^). 

The ratio /m = 65 is one expression of DCSB in the parametrisation of S{p). It 
emphasises the dramatic enhancement of the dressed-quark mass function at infrared momenta. 
Another is the chiral-limit vacuum quark condensate [11] 



(5-'?)? = A^T^^ln-^, (A.25) 



which assumes the value (/Iqcd = 0.2 GeV) 

- (w)?=i Gcv = (0.221 GeV)^ (A.26) 
A detailed discussion of the vacuum quark condensate in QCD can be found in Ref. [701 [71] 



A. 2. 2 Diquark Bethe-Salpeter amplitudes 

The rainbow-ladder DSE truncation yields asymptotic diquark states in the strong interaction 
spectrum. Such states are not observed and their appearance is an artefact of the truncation. 
Higher-order terms in the quark-quark scattering kernel, whose analogue in the quark-antiquark 
channel do not much affect the properties of vector and flavour non-singlet pseudoscalar mesons, 
ensure that QCD's quark-quark scattering matrix does not exhibit singularities which corre- 
spond to asymptotic diquark states [TH]. Nevertheless, studies with kernels that don't generate 
diquark bound states, do support a physical interpretation of the masses, 'm(qq)jp, obtained 
using the rainbow-ladder truncation: the quantity l(qq)jp ~ l/m(gg)^^ may be interpreted as 
a range over which the diquark correlation can propagate within a baryon. These observations 
motivate an Ansatz for the quark-quark scattering matrix that is employed in deriving the 
Faddeev equation: 

[Mqq{k,q-K)]l.:= Yl r'\k--K)A-'\K)r'\q-K). (A.27) 

JP=0+ ,1+ 

= 10"'' in Eq. (|A.19[) acts only to decouple the large- and intermediate-p^ domains. 
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One manner of specifying the in Eq. (|A.27[l is to employ the solutions of a rainbow- 
ladder quark-quark Bethe-Salpeter equation (BSE), as e.g. in Refs. [231 1271 138) . Using the 

properties of the Gell-Mann matrices one finds easily that Fq := F C' satisfies exactly the 
same equation as the colour-singlet meson but for a halving of the coupling strength [25) . 
This makes clear that the interaction in the 3c (qq) channel is strong and attractive. 

A solution of the BSE equation requires a simultaneous solution of the quark-DSE. How- 
ever, since we choose to simplify the calculations by parametrising S{p), we also employ that 
expedient with F'^ , using the following one-parameter forms: 

F°*{k;K) = -l-H^Ci'ysir2HkV^o+), (A.28) 
t'r^^^(fc;if) = -l^H^i^^Ct':F{k^/ul+), (A.29) 

with the normalisation, A/""^ , fixed by requiring 

(A.30) 

d^q 



2K^ 



d 



/i4 
^ r(q; ~K) S{q + Q/2) F{q; K) S'^{~q + Q/2). (A.31) 

The Ansatze of Eqs. ()A.28|I . ()A.29)) retain only that single Dirac-amplitude which would 
represent a point particle with the given quantum numbers in a local Lagrangian density. They 
are usually the dominant amplitudes in a solution of the rainbow-ladder BSE for the lowest 
mass diquarks [26l[27] and mesons p0l[72l[73] . 



A. 2. 3 Diquark propagators 

Solving for the quark-quark scattering matrix using the rainbow-ladder truncation yields free 
particle propagators for A'-' in Eq. ()A.27|I . As already noted, however, higher-order contri- 
butions remedy that defect, eliminating asymptotic diquark states from the spectrum. The 
attendant modification of A'^ can be modelled efficiently by simple functions that are free- 
particle- like at spacelike momenta but pole- free on the timelike axis [18]; namelylf] 

A''\k) = -^:F{K'/u;1+), (A.32) 

KUk) = (s,^ + ^^) ^TiK'/-^!^), (A.33) 

where the two parameters rrijp are diquark pseudoparticle masses and lJjp are widths charac- 
terising F'^ . Herein we require additionally that 



bjjp — i rn^jp , (A. 34) 



which is a normalisation that accentuates the free-particle-like propagation characteristics of 
the diquarks within the hadron. 

Appendix B: Euclidean Conventions 

In our Euclidean formulation: 

4 

p-q = '^Piqi; (B.i) 

1=1 

{7m.7^} = 2(5p^ ; 7i = 7M; o-m-- = ^ [7m> 7>'] ; tr [757^7^ 7p7<t] = -4 e^^p„ , £1234 = 1 . (B.2) 



^These forms satisfy a sufficient condition for confinement because of the associated violation 
of reflection positivity. See Sect. 2 of Ref. [14] for a brief discussion. 
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A positive energy spinor satisfies 

u{P, s){i'y ■ P + M) = = {iy ■ P + M) u{P, s) 
where s = ± is the spin label. It is normalised: 

u{P, s) u{P, s) = 2M , 

and may be expressed explicitly: 



i(P, s) = VAf - i£ 



with £ = WP'^ + 



rr P 



For the free-particle spinor, u{P, s) = u(P, 3)^74. 

The spinor can be used to construct a positive energy projection operator: 



A, 



^= ^ E "(^' ^(^' ^) = ^ (-^^ -P + M) 



2M 
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(B.3) 



(B.4) 



(B.5) 



(B.6) 



(B.7) 



A negative energy spinor satisfies 

v{P, s) (i-f-P- M) = = (i7 • P - M) v{P, s) , 



(B.i 



and possesses properties and satisfies constraints obtained via obvious analogy with u{P, s) 
A charge-conjugated Bethe-Salpeter amplitude is obtained via 



r{k-P) =C'' r{-k;Pf C , 



(B.9) 



where "T" denotes a transposing of all matrix indices and C — 7274 is the charge conjugation 
matrix, = — C. 



Appendix C: Nucleon-Photon Vertex 

In order to explicate the vertex depicted in Fig. [2] we write the scalar and axial- vector compo- 
nents of the nucleons' Faddeev amplitudes in the form [cf. Eq. (IA.11|) ] 



For explicit calculations, we work in the Breit frame: P^ — P^^ — Qn/2, P^ = P^f ^ + Qft/2 
and P^^ — (0, 0, 0, i^/mI + QV^), and write the electromagnetic current matrix element as 
[cf. Eq. ©] 



P'Uf. 



A+{P') 



'y^iGE + Mn-^ — [Ge — Gm) 



p2 

^BF 



A+{P), 



Wf W)i ^"^'^-^ ^) ■ 



(C.2) 
(C.3) 



In Fig. [2] we have separated the current, Jf^ {p, P' \k, P), into a sum of six terms, each of 
which we subsequently make precise. NB. Diagrams 1, 2 and 4 are one-loop integrals, which we 
evaluate by Gaui3ian quadrature. The remainder. Diagrams 3, 5 and 6, are two-loop integrals, 
for whose evaluation Monte-Carlo methods are employed. A technical aspect concerning the 
computation is described in [Appendix D:| Chebyshev expansion. 
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C.l Diagram 1 

This represents the photon couphng directly to the bystander quark. It is expressed as 

jr = S(p,)fr(p,; k,)S{k,) (z\°^ {ks) + A'* (ks)) {2^f5\p -k-f,Q), (C.4) 

where /^"(p,; kq) = Qq r^{pq; kq), with Qq = diag[2/3, —1/3] being the quark electric charge 
matrix, and r^{pq; kq) is the dressed-quark-photon vertex. In Eq. (|C.4|) the momenta are 

kq = rjP + k, Pq = rjP' +p, 

kd ^ fjP ~ k , pd = fjP' - p, 

with rj + fi = 1. The results reported herein were obtained with rj = 1/3, which provides a single 
quark with one-third of the baryon's total momentum and is thus a natural choice. Notably, 
as our approach is manifestly Poincare covariant, the precise value is immaterial so long as 
the numerical methods preserve that covariance. Calculations converge most quickly with the 
natural choice. 

It is a necessary condition for current conservation that the quark-photon vertex satisfy 
the Ward-Takahashi identity: 

Q^ir^{ei,e2) = s-'iii) ^ s-\£2) , (c.6) 

where Q = — ^2 is the photon momentum flowing into the vertex. Since the quark is dressed. 
Sec. I A. 2. n the vertex is not bare; i.e., 1)1(^1,^2) 7^ 7/j. It can be obtained by solving an inho- 
mogeneous Bethe-Salpeter equation, which was the procedure adopted in the DSE calculation 
that successfully predicted the electromagnetic pion form factor |39l I73j . However, since we 
have parametrised S{p), we follow Ref. [11] and write [74] 

ir^{ei,£2) = iSA{fi,el) tm + 2fc^ [n • k^ Aa{i\,iI) + AB{i\,il% (c.7) 

with = (^1 ^2)/2, Q = {ti- £2) and 



SF{el£l)^^[Fiel)+F{£l)], = (C.8) 

where F = A,B; viz., the scalar functions in Eq. HA.17|) . It is critical that J), in Eq. (|C.7|) satisfies 
Eq. (|C.6|) and very useful that it is completely determined by the dressed-quark propagator. 



C.2 Diagram 2 

This figure depicts the photon coupling directly to a diquark correlation. It is expressed as 

jdq 



J,7 = A 



\Pd) [r^''{pd;kd)Y A^{kd)S{kq){2nf5\p~k + rjQ), (C.9) 

with [f^'ipd-, kd)Y' = diag[Qo+rO+, Qi+r^i^], where Qo+ = 1/3 and pf is given in Eq. (ICT4ll . 
and (5i+ = diag[4/3, 1/3, —2/3] with P^ given in Eq. (|C.16|) . Naturally, the diquark propaga- 
tors match the line to which they are attached. 

In the case of a scalar correlation, the general form of the diquark-photon vertex is 

{£i,e2) = 2k^f+[k\k- Q, q") + /_ (fc^ fc • Q, q") . (c.io) 

If one is dealing with an elementary scalar correlation, then the Ward-Takahashi identity reads: 

Q^p°^i£i,£2) = n°^i£l)-n°*{£l), n-'"" {£^) = {A-''' {£^)}-\ (C.ll) 

However, for a composite system of the type we consider this identity is modified; viz. [75], 



hF, 



fi£i,£2) = [n'>\£l) - n"U£l)] FqqiQ'), (C.12) 



where 



1 + e'^i^Q 
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is a form factor describing the distribution of charge within the correlation. 

The evaluation of scalar diquark elastic electromagnetic form factors in Ref. [38] is a first 
step toward calculating (^1,^2)- However, in providing only an on-shell component, it is 
insufficient for our requirements. We choose to adapt Eq. ()C.7|) to our needs and employ 



(4, ^2) = fcp ^no+ i^ijl) F,,{Q^) , (C.14) 

with the definition of '^^0+ (^i) ^2) apparent from Eq. HC.8|I and the value of rqq given in Eq. (|10|) . 

Equation (|C.14|l is an Ansatz that satisfies Eq. (|C.12[l . is completely determined by quan- 
tities introduced already and is free of kinematic singularities on the relevant domain. It im- 
plements /_ = 0, which is a requirement for elastic form factors, and guarantees a valid nor- 
malisation of electric charge; viz.. 

Urn rf {e',i) = 21,^ n"^ {f) 2 i, , (c.15) 

owing to Eq. (IA.34p . NB. We have factored the fractional diquark charge, which therefore ap- 
pears subsequently in our calculations as a simple multiplicative factor. 

For the case in which the struck diquark correlation is axial-vector and the scattering is 
elastic, the vertex assumes the form [76) 

3 
i = l 

with {Tc,fi{£) = Sc,f,-£o.£f3/£'') 

r!!U£i,e2) = {h+l2)^T^x{h)Tx0{l2)F^{tl,ll), (C.17) 

^M'i/3(^1'^2) = [T^o.{h)T,3p{l2)tlp+T^p{i2)T^p{il)i2p]F2(A,il), (CIS) 
rflp{il,i2) = -^-^{h+i2)^.T^p{h)l2pTpx{l2)hxF3{l\,ll). (C.19) 

This vertex satisfies: 

tic Fl^^f, (^1 , ^2 ) = = r,l^a {h , £2 ) £213 , (C.20) 

which is a general requirement of the elastic electromagnetic vertex of axial- vector bound states 
and guarantees that the interaction does not induce a pseudoscalar component in the axial- 
vector correlation. We note that the electric, magnetic and quadrupole form factors of an 
ax;ial- vector bound state are expressed [76] 

(Q^) = Fi + I G'q (Q^) , n+ = , (C.21) 

G'm{Q^) ^ -F2{Q^) , (C.22) 
Gq'(O') =Fi(Q')+F2(Q') + (l + ri+)F3(Q'). (C.23) 

Owing to the fact that Fq :— F' C satisfies exactly the same Bethe-Salpeter equation as 
the colour-singlet meson but for a halving of the coupling strength, the vector meson form 
factor calculation in Ref. [40] might become useful as a guide in understanding the form factors 
in Eqs. HC.16|I - (|C.19|I . However, in providing only an on-shell component, that information is 
insufficient for our requirements. Hence we employ the following Ansatze: 

Fr{i\,ll) = A^,+ (£lll)Fqq{Q^), (C.24) 
F2{£l,£l) = -Fi-h(l-ri+)(ri+Fi + l-^i+)d(ri+), (C.25) 
F3{£l£l) = -(xi+(l-ri+)d(ri+)+Fi+F2)rf(ri+), (C.26) 



''If the scattering is inelastic the general form of the vertex involves eight scalar functions [77] . 
For simplicity, we ignore the additional structure in this Ansatz. 
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with d(x) — 1/(1 + x)^. This construction ensures a vahd electric charge normalisation for the 
axial- vector correlation; viz., 

Urn r'^tfsie'J) = T^p{e) n'* (f) T.^f) , (c.27) 

owing to Eq. (|X34)) . and current conservation 

lim Qf.r^lB{ei,i2)^0. (C.28) 
The diquark's static electromagnetic properties follow: 

(0) = 1 , G'm (0) = Mi+ , G'q (0) = -xi+ . (C.29) 

For an on-shell or pointlike axial-vector: = 2; and Xi+ ~ 1- In addition, Eqs. HC.16|I - (|C.19|) 
with Eqs. (|C.24[| - (|C.26|I realise the constraints of Ref. [7S]; namely, independent of the values 
of fj.i+ & Xi+ 1 the form factors assume the ratios 

GIUq') ■■ G'm{Q') ■■ G'q{Q') '^'=°° (1 - |ri+) : 2 : -1. (C.30) 

It is noteworthy that within a nucleon the diquark correlation is not on-shell. Hence, in 
contrast with Ref. [31], we do not assume herein that a point-particle value for the magnetic 
moment in Eq. (|C.29|) serves as a good reference point. Instead we employ the value determined 
in Ref. [33]: 

= 0.37, (C.31) 

which is in accord with that obtained following the approach in Ref. [23]- While equally one 
need not employ the point-particle value for Xi+ i changing to Xi+ = has little impact on the 
results [31]. We therefore stay with Xi+ ~ 1- 



C.3 Diagram 3 

This image depicts a photon coupling to the quark that is exchanged as one diquark breaks up 
and another is formed. It is expressed as 

= -ls{k,)A\ka)r\p^,ka)S^iq)rr^iq',q)S^iq')r'^ip'2,p,)A^p,)S{p,) , (C.32) 

wherein the vertex appeared in Eq. (|C.4[l . While this is the first two- loop diagram we have 
described, no new elements appear in its specification: the dressed-quark-photon vertex was 
discussed in Sec. IC.ll In Eq. HC.32P the momenta are 

q = rjP — 77P' — p ~ k , q' = vP' — rjP — p ~ k , 

Pi = {p,-q)/2, p'2 = (-k, + q')/2, (C.33) 

p[ = {pg-q')/2, p2^{-k, + q)/2. 

It is noteworthy that the process of quark exchange provides the attraction necessary in 
the Faddeev equation to bind the nucleon. It also guarantees that the Faddeev amplitude has 
the correct antisymmetry under the exchange of any two dressed-quarks. This key feature is 
absent in models with elementary (noncomposite) diquarks. The full contribution is obtained 
by summing over the superscripts which can each take the values 0^, 1^. 



C.4 Diagram 4 

This differs from Diagram 2 in expressing the contribution to the nucleons' form factors owing 
to an electromagnetically induced transition between scalar and axial- vector diquarks. The 
explicit expression is given by Eq. HC.9|) with [r^'^{pd\ fcd)]'"^ = 0, and [r^''{pd\ fcd)]^'^ = Psa 
and [i7i'(Pdi ^d)]^'^ = Pas- This transition vertex is a rank-2 pseudotensor, kindred to the 
matrix element describing the p^'ix^ transition [79], and can therefore be expressed 

rj^(^l,^2) = -PTs{ilA2) = ^T(£i,^2)e^.pAW2A, (C.34) 
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where 7, a. are, respectively, the vector indices of the photon and axial-vector diquark. For 
simplicity we proceed under the assumption that 

T{h,h)=t^r\ (C.35) 

A typical on-shell value for the dimensionless normalisation is Kr ~ 2 [SD]. However, as with 
, we recognise herein that this value is not a useful reference point because, for the processes 
described by Fig. (2] nr can be much smaller in magnitude. We use the value determined in 
Ref. [33]: 

KT = 0.12. (C.36) 
This diagram impacts upon the nucleons' magnetic form factors pi 1311 155] . 

C.5 Diagrams 5 &: 6 

These two-loop diagrams are the so-called "seagull" terms, which appear as partners to Dia- 
gram 3 and arise because binding in the nucleons' Faddeev equations is effected by the exchange 
of a dressed-quark between nonpointlike diquark correlations [36] . The explicit expression for 
their contribution to the nucleons' form factors is 

- r{pi,k,)S^{q)Xf,{-k„-q,pa)) A'{pd)S{p,), (C.37) 

where, again, the superscripts are summed. 

The new elements in these diagrams are the couplings of a photon to two dressed-quarks as 
they either separate from (Diagram 5) or combine to form (Diagram 6) a diquark correlation. As 
such they are components of the five point Schwinger function which describes the coupling of a 
photon to the quark-quark scattering kernel. This Schwinger function could be calculated, as is 
evident from the computation of analogous Schwinger functions relevant to meson observables 
[81] . However, such a calculation provides relevant input only when a uniform truncation of the 
DSEs has been employed to calculate each of the elements described hitherto. We must instead 
employ an algebraic parametrisation [36], which for Diagram 5 reads 

xi"{k^Q) = [r'\k ^ Q/2) - r'\k)\ 

with k the relative momentum between the quarks in the initial diquark, eby the electric charge 
of the quark which becomes the bystander, and ecx the charge of the quark that is reabsorbed 
into the final diquark. Diagram 6 has 

xf{k,Q) = e^,^^±f^^[r'\k + Q/2)-r'\k)] 

+ -4S^^[^-^'(^-«/2)-^''(^^)]' (^-^9) 

where t'' {t} is the charge-conjugated amplitude, Eq. (]B.9|I . Plainly, these terms vanish if the 
diquark correlation is represented by a momentum-independent Bethe-Salpeter-like amplitude; 
i.e., the diquark is pointlike. 

It is naturally possible to use more complicated Ansatze [23]. However, like Eq. ()C.14[) . 
Eqs. (]C.38[) & (]C.39P are simple forms, free of kinematic singularities and sufficient to ensure 
the nucleon-photon vertex satisfies the Ward-Takahashi identity when the composite nucleon 
is obtained from the Faddeev equation. 

Appendix D: Chebyshev expansion 

In solving the Faddeev equation we employ a Chebyshev expansion of the scalar functions 
appearing in the Faddeev amplitude and wave function in order to restrain the use of computer 
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memory. (See, e.g., Ref. [20].) The results herein were obtained with twelve terms in both. The 
Chebyshev-expanded functions then define the Faddeev amplitude that appears and is evaluated 
in the expressions for the form factors. Without due care, this can lead to a problem; namely, 
with increasing a function can be computed outside the expansion's domain of convergence. 

Consider a function F{k^, k-P;P^), which represents a term in the Faddeev amplitude. It is 
a function of only two variables: k^ and k ■ P, where k is the relative quark-diquark momentum, 
because the total momentum always satisfies P^ = —M^, where M is the bound-state's mass. 
In the bound-state's rest frame one can define an angle a through 

zlfclA/cosQ — A: • P. (D.l) 

Then, with {Ui{x), j = 1 ... oo} being Chebyshev polynomials of the second kind, 

F{k'^,k- P;-M'^) ^ lim V ■'F(|fcj,a/; -A/^) C/j(cos/3) . (D.2) 

Nm^OO ^ — ^ 
3 = 

For any finite Nm the expansion in Eq. (|D.2|I is a true approximation to the k ■ P-dependence 
of the function F in the sense that, with increasing Nm, the right-hand-side (rhs) is uni- 
formly pointwise an increasingly accurate representation of the function. The Ihs of Eq. (|D.2|I 
is Poincare invariant. Hence, in the limit Nm — > oo, so is the rhs. These statements are true so 
long as cos Of defined in Eq. (|D.1|) satisfies — 1 < cosq < 1. 

In calculating a form factor one must compute the Faddeev amplitude of a bound-state 
that is not at rest. In the Breit frame, e.g., the total momentum can be written as P = 
(0,0, ±Q/2,iF((3/2)), where E^{Q/2) = A/^ -I- Q^/4, the bound-state is moving with three 
momentum ±Q/2 and 

k- P = ±^\k\Q cose sin (3 + i\k\E{Q) cos (3 , (D.3) 

with k expressed using the standard definition of hyperspherical coordinates. In principle, as 
demonstrated in Ref. 00], this is not a problem in a Poincare covariant framework. However, it 
can consume large amounts of computer memory and time. We therefore proceed by writing 

k - P = i\k\E{Q) 

in which case the real and imaginary parts of z are bounded in magnitude by one, and then 
define 

P(fc^ fc • (P ± Q/2); -M^) = J2 'F{\k\,iE{Q)--M^) U,{z) . (D.5) 
Appendix E: Form Factor Notation 

We represent all form factors by their usual symbols. Hence, the notation can be introduced 
via an exemplar; viz., the proton's Pauli form factor, Ff . 

• Pf - Sum of all contributions to Ff that can be represented by Diagram 1 in Fig. [2] i.e., 
in which the photon interacts with a bystander quark, either u or d. Pf = F['''{Q^ — 0) 
gauges the probability that the photon interacts with a bystander quark. 

• P^'" - Sum of all contributions to Pf that can be represented by either Diagram 2 or 
4; i.e., in which the photon interacts with a diquark correlation, either scalar or axial- 
vector, or excites a transition between them. Pf — Ff'''{Q^ = 0) gauges the probability 
that the photon interacts with a diquark. 

• Ff''^ - Sum of all contributions to Pf that can be represented by one of Diagrams 3, 5 
or 6; i.e., in which the photon interacts with a diquark in association with its breakup. 
pp,e _ pP'''(^Q^ — 0) gauges the probability that the photon acts in association with 
diquark breakup. 

NB. Pf" + Pf + Pf = Pf . 

• p^'" - Sum of all contributions to P^ in Fig. [2] that are proportional to the charge of a 
u-quark, Eu; i.e., the total li-quark contribution F[ . 



iQ 



2E{Q) 



cos e sin P + cos (3 



i\k\E{Q)z, 



(D.4) 
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• F^''^'^ ~ Sum of all contributions to Ff'^ that can be represented by Diagram 1 in Fig. [J] 
i.e., in which the photon interacts with a bystander u-quark. 

• Ff''^'^ - Sum of all contributions to Ff that can be represented by either Diagram 2 or 4 
and are proportional to e^, i.e., in which the photon resolves a it-quark within a diquark 
correlation. 

• Ff''^'^ - Sum of all contributions to _Ff that can be represented by one of Diagrams 3, 
5 or 6 and are proportional to e„; i.e., in which the photon interacts with a u-quark in 
association with the breakup of a diquark. 

NB. i^f'"'" + Ff'"'" + Ff'"'" = Ff'"; Ff'"(0) = 2e„; 2e„Pf'"'" := Ff'°'"(Q2 ^ o), 
a — q,d, e. 

• Ff ''^ and related functions are defined in direct analogy with those connected to Ff'". 
NB. Ff^''''' + Ff'''-'' + Ff = Ff'''; Ff''*(0) = e^; eaPt"''' - F^^-^Q^ = 0), a = g,d,e. 

• Ff '° - Sum of all contributions to Ff in Fig. [2] that involve a scalar diquark component 
in both 'Fi and li//. Pf'" = Ff'"{Q^ = 0) gauges the probability that the photon interacts 
with a scalar diquark component of the nucleon. 

• Ff - Sum of all contributions to Ff that involve an axial-vector diquark component in 
both 'Pi and 'Pf. Pf''^ = Fl''^{Q^ = 0) gauges the probability that the photon interacts 
with an axial- vector diquark component of the nucleon. 

• Ff'"^ - Sum of all contributions to Ff in which the diquark component of 'Fi is different 
to that in 'Pf. P^'"" = F^'"^{Q = 0) gauges the probability that the photon induces a 
transition between diquark components of the incoming and outgoing nucleon. 

NB. Ff'"-hFf '" 4- Ff ^-^Ff. 

• Ff - Sum of all contributions to F[ in Fig.[2]that involve a scalar diquark component 
in both 'Pi and 'Pf, and are proportional to e^; i.e., in which a w-quark is resolved in the 
presence of a scalar diquark. 

• Ff - Sum of all contributions to Ff that involve an axial- vector diquark component 
in both 'Pi and 'Pf, and are proportional to e„. 

• Ff''"'^ - Sum of all contributions to Ff that are proportional to Cu and in which the 
diquark component of 'Pi is different to that in 'Pf. 

NB. Ff'"'" + Ff'"'" + Ff"'" = Ff "; 2e„Pf"'" — Ff'"'"(Q2 = 0), a = s,a,m. 

• Ff and similar functions are defined in direct analogy with those connected to Ff 
NB. Ff + Ff ''^ + Ff = Ff*; e^Pf - Ff '"■''(Q^ = Q), a = s, a, m. 
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